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Abstract 

Let IK be an ordinary differential field with derivation d. Let V he a, 
system of n linear differential polynomial parametric equations in n — 1 
differential parameters with implicit ideal ID. Given a nonzero linear 
differential polynomial A in ID we give necessary and sufficient condi- 
tions on A for 'P to be n — 1 dimensional. We prove the existence of 
a linear perturbation of V so that the linear complete differential 
resultant 9CRes0 associated to is nonzero. A nonzero linear dif- 
ferential polynomial in ID is obtained from the lowest degree term of 
dCRes^ and used to provide an implicitization algorithm for V. 



1 Introduction 

The use of algebraic elimination techniques such as Groebner basis and mul- 
tivariate resultants to obtain the implicit equation of a unirational algebraic 
variety is well known (see for instance [4], [5]). The development of similar 
techniques in the differential case is an active field of research. In [7] character- 
istic set methods were used to solve the differential implicitization problem for 
differential rational parametric equations and alternative methods are emerg- 
ing to treat special cases. In [13], linear complete differential resultants were 
used to compute the implicit equation of a set of linear differential polynomial 
parametric equations (linear DPPEs). As in the algebraic case differential re- 
sultants often vanish under specialization and we are left with no candidate 
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to be the implicit equation. This reason prevented us from giving an algo- 
rithm for differential implicitization in [13]. Motivated by Canny's method in 
[2] and its generalizations in [6] and [10] in the present work we consider a 
linear perturbation of a given system of linear DPPEs and use linear complete 
differential resultants to give a candidate to be the implicit equation of the 
system. 

Given a system V{X, U) of n linear ordinary differential polynomial para- 
metric equations Xi = Pi{U), . . . , a:„ = Pn{U) in — 1 differential parameters 
Ml, ... , Un-i (we give a precise statement of the problem in Section 2) we give 
an algorithm to decide if the dimension of the implicit ideal ID of P is n — 1 
and in the affirmative case provide the implicit equation of V. 

The linear complete differential resultant dCRes{xi—Pi{U) , . . . , Xn—PniU)) 
is the algebraic resultant of Macaulay of a set of differential polynomials with 
L elements. It was defined in [13] as a generalization of Carra'Ferro's dif- 
ferential resultant [3] (in the linear case) in order to adjust the number L of 
differential polynomials to the order of derivation of the variables Mi, . . . , m„ in 
Fi = Xi-P^{U). 

In this paper, we provide a perturbation '^^(X, U) of P(X, U) so that the 
linear differential polynomials Fi — p0i(f/), . . . , Fn—p4>n{U) have nonzero lin- 
ear complete differential resultant dCIies^{p) which is a polynomial depend- 
ing on p. It will be shown that the coefficient of the lowest degree term of 
(9CRes0(p) is a nonzero linear differential polynomial that belongs to the im- 
phcit ideal ID of V{X,U). In fact, if 9CRes0(p) has a nonzero independent 
term it equals c)CRes(Fi, . . . , F„) and as proved in [13] it gives the implicit 
equation of V{X, U). 

The main result of this paper generalizes the result previously mentioned 
from [13]. Given a nonzero linear differential polynomial A in ID necessary 
and sufficient conditions on A are provided so that A{X) = is the implicit 
equation of V{X,U). The higher order terms in the equations of V{X,U) 
and the rank of the coefficient matrix of the set of L polynomials used to 
construct the differential resultant c}CRes(Fi, . . . , F„) play a significant role in 
this theory. The fact that we are dealing with linear differential polynomials 
will be also relevant, allowing us to treat then using differential operators. 

The paper is organized as follows. In Section 2 we introduce the main 
notions and notation. Next we review the definition of the linear complete dif- 
ferential (homogeneous) resultant in Section 3. The main definitions regarding 
linear differential polynomials in ID are given in Section 4. The next section 
contains the main result of the paper, namely a characterization of ID in the 
n—1 dimensional case is provided in Section 5. In Section 6 we give a pertur- 
bation V^{X, U) of V{X, U) with nonzero differential resultant and use it to 



2 



obtain a nonzero linear differential polynomial in ID candidate to provide the 
implicit equation. Finally in Section 7 we give the implicitization algorithm 
and examples. 

2 Basic notions and notation 

In this section, we introduce the basic notions related to the problem we deal 
with (as in [12] and [13]), as well as notation and terminology used throughout 
the paper. For further concepts and results on differential algebra we refer to 
[8] and [9]. 

Let K be an ordinary differential field with derivation 9, ( e.g. Q(t), d = 
^). Let X = {xi,...,Xn} and U = . . . , be sets of differential 

indeterminates over K. Let No = {0,1,2,...,?t,,...}. For A; G No we denote 
by Xik the k-th derivative of Xj, for x^o we simply write Zj. Given a set Y of 
differential indeterminates over IK we denote by {Y} the set of derivatives of 
the elements of Y, {Y} = {d^y \ y E Y, /c G No}, and by ]K{X} the ring of 
differential polynomials in the differential indeterminates xi, . . . that is 

= K[xik I z = l,...,n, ke No]. 

Analogously for Kjf/}. 

As defined in [13] we consider the system of linear DPPEs 

r xi = Pi(f/) 

V{X,U) = \ : (1) 

[ Xn = Pn{U). 

where Pi,...,Pn G K{f/} with degree at most 1 and not all Pj G K, z = 
l,...,n. There exists G IK and an homogeneous differential polynomial 
Hi G K{U} such that 

F,(X, U)=x,- P,{U) =x,-a, + H,{U). 

Given P G K{X UU} and ?/ G X U f/, we denote by ord(P, y) the order of P 
in the variable y. If P does not have a term in y then we define ord(P, y) = —1. 
To ensure that the number of parameters is n — 1, we assume that for each 
jG{l,...,n — 1} there exists i G {1, . . . ,n} such that ord(Pi, Uj) > 0. 

The implicit ideal of the system (1) is the differential prime ideal 

ID = {/ G K{X} I /(Pi(f/), . . . , PniU)) = 0}. 

Given a characteristic set C of ID then n— | C | is the (differential) dimension 
of ID, by abuse of notation, we will also speak about the dimension of a DPPEs 
system meaning the dimension of its implicit ideal. 
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If dim(ID) = n — 1, then C = {A{X)} for some irreducible differential 
polynomial A G K{X}. The polynomial A is called a characteristic polynomial 
of ID. The implicit equation of a (ra — l)-dimensional system of DPPEs, in 
n differential indeterminates X = {xi, . . . , x„}, is defined as the equation 
A{X) = 0, where A is any characteristic polynomial of the implicit ideal ID of 
the system. 

Let K.[d] be the ring of differential operators with coefficients in K. If K is 
not a field of constants with respect to d, then ]K[d] is not commutative but 
dk — kd = d{k) for all /c G K. The ring K.[d] of differential operators with 
coefficients in K is left euclidean (and also right euclidean). Given C, C G IK[9], 
by applying the left division algorithm we obtain q,r & ^[d], the left quotient 
and the left reminder of £ and C respectively, such that £ = £'g + r where 
deg(r) < deg(£'). 

3 Linear complete differential resultants 

We review next the results on linear complete differential resultants from [13] 
that will be used in this paper. 

Let © be a differential integral domain, and let G D{f/} be a linear ordi- 
nary differential polynomial of order Oi, i = 1, . . . ,n. For each j G {1, . . . , n— 1} 
we define the positive integers 



Let = then the completeness index 7(/i, . . . , fn) verifies 7(/i, . . . , fn) < 

— Oi, for alH G {1, . . . , n}. 

The linear complete differential resultant c?CRes(/i, . . . , fn) is the Macaulay's 
algebraic resultant of the differential polynomial set 



The set PS(/i, . . . , /„) contains L = J2^=ii^ — Oj — 7 + 1) polynomials in the 
following set V of L — 1 differential variables 

V = {uj,Uji . . .^UjN-^j-y I Ij = Ijifi, ■ ■ ■,fn),j = l,...,n-l}. 

Let hi G D{f/} be a linear ordinary differential homogeneous polynomial 
of order Oj, i = 1, . . . ,n with = Y17=i > 1- We define the differential 



7i(/i: ■■■Jn) ■■= min{oi - ord(/i,Mj) | i G {1, . . . 



n-l 




PS(/i,...,/„) : = 

. . . , 9/„ /, h = 1, . . . , n, 7 = 7(/i, . . . , /n)}. 
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polynomial set 

PS\h, ...,hn):= . . . , dh„ hi I 

ie{l,...,n}, iV-Oi-7-l>0, 7 = 7(/ii,...,/i„)}. 

Observe that > 1 implies PS^ {hi, hn) 7^ 0. The linear complete dif- 
ferential homogenous resultant 9CRes'*(/ii, . . . , hn) is the Macaulay's algebraic 
resultant of the set PS'^(/ii, . . . , hn)- The set PS'^(/ii, . . . , hn) contains = 
Sr=i(^ ^ Oj — 7) polynomials in the set V'* of L'^ differential variables 

V'^ = {uj,Uji . . . ,UjN-^j-i-i I 7j = Ijihi, hn),j = 1, . . . ,n - 1}. 

We review next the matrices that will allow the use of determinants to 
compute dCRes{ fi, fn) and dCRes^ {hi, hn)- The order ui < ■■■ < 
induces an orderly ranking on U (i.e. and order on {U}) as follows (see 
[8], page 75): u < du, u < u* ^ du < du* and k < k* ^ d^u < d^* u* , for all 
n, u* e U, k, k* e Nq. We set 1 < Ui. 

For i = 1, . . . , n, 7 = 7(/i, . . . , /„) and k = 0, . . . , — Oj — 7 define de 
positive integer l{i, k) = {i — 1){N — 7) — J2^h=i Oi + i + k in {1, ... , L}. The 
complete difFerential resultant matrix M{L) is the L x L matrix containing the 
coefficients of d^~°^~"'~^ fi as a polynomial in D[V] in the l{i, k)-i\i row, where 
the coefficients are written in decreasing order with respect to the orderly 
ranking on U . In this situation: 

9CRes(/i,...,/„) = det(M(L)). 

If iV > 1, for 7 = 7(/ii, . . . , hn), i e {I, . . . ,n} such that iV-Oi -7-I >0 
and k = 0, . . . , — Oj — 7 — 1 define the positive integer l^{i, k) = {i — 
l)(iV - 7 - 1) - X;l=\ Oi + i + k in {!,..., L^}. The complete differential 
homogeneous resultant matrix M{L'^) is the L'^ x matrix containing the 
coefficients of d^~°^~"'~^~^hi as a polynomial in D[V''] in the l^{i,k)-t\i row, 
where the coefficients are written in decreasing order with respect to the orderly 
ranking on U . In this situation: 

aCRes'^(/ii, ...,hn) = det{M{L^)). 

Throughout the remaining parts of the paper we will say differential ( homo- 
geneous) resultant always meaning linear complete differential ( homogeneous) 
resultant. 

3.1 Linear complete differential resultants from linear 
DPPEs 

We highlight in this section some facts on differential resultants of the differ- 
ential polynomials Fi and Hi obtained from a system of linear DPPEs as in 
Section 2. 
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Let 7 = 7(Fi, . . . , F„) = 7(^/1, . . . , and © = K{X}. Then aCRes(Fi, . . . , 
and dCRes^{Hi, . . . , are closely related as shown in [13], Section 5. Since 
Fi{X,U) = Xi - ai + Hi{U), if iV > 1 the matrix M{L'^) is a submatrix of 
M[L) obtained removing n specific rows and columns. This fact together with 
the identities below allowed to prove that (when > 1) 

aCRes(Fi, . . . , F„) = ^ dCRes^{Hi, . . . , = 0. 

The next matrices will play an important role in the remaining parts of the 
paper. 

• Let S be the n x (n — 1) matrix whose entry (2, j) is the coefficient of 
Un-jo.-'Yj in Fi, i G {1, . . . ,n}, j G {1, . . . , n — 1}. We call S the leading 
matrix of V{X, U). For i G {1, . . . n} let Si be the (n — 1) x (n — 1) matrix 
obtained by removing the i-th row of S. 

• Let Mi_i be the L x (L — 1) principal submatrix of M (L). We call M^-i 
the principal matrix of V{X, U). 

Let X = {xi, Xii, . . . , XiN-o,-'y M = 1, • • • 5 n}. Given x E X, say x = Xik 
with /c G {0, 1, . . . , A^— Oj— 7}, let us call the submatrix of Ml-i obtained by 
removing the row corresponding to the coefficients of d'^Fi = xik + d^{Hi{U) — 
Qi). Then, developing the determinant of M{L) by the last column we obtain 

n N—Oi—'y 

aCRes(Fi,...,F„) = ^ hkdei{M^J{x,k-d''ai), (2) 

with bik = ±1 according to the row index of Xik — d^ai in the matrix M{L). 
Also, for every i G {1, . . . , n}, there exists a G N such that 

det(M,,^_„^_^) = (-l)»aCRes^(i/i, . . . , det(S',). (3) 

4 The implicit ideal ID 

Let V{X, U), Fi, Hi be as in Section 2. Let PS = PS(Fi, . . . , F„) and let ID 
be the implicit ideal of V{X,U). In this section, we review the computation 
of ID in terms of characteristic sets (see [7] and [13]) and give some definitions 
related with linear differential polynomials in ID that will be important in the 
remaining parts of the paper. 

Let X and V be as in Section 3.1 and observe that PS C K[A'][V]. Let (PS) 
be the ideal generated by PS in K[X] [V] and let [PS] be the differential ideal 
generated by PS in K{X}. 

The order Xn < ■ ■ ■ < Xi induces a ranking on X as follows (see [8], page 
75): X < dx and x < x* ^ d^x < d^* x* , for all x, x* G X, fc, k* G Nq. 
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• We call TZ the ranking on X UU that eliminates X with respect to U, 
that is d'^x > d''*u, for all x E X, u E U and k, k* E Nq. 

• We call TZ* the ranking on X U [/ that eliminates U with respect to X, 
that is d^x < d^* u, for al\ x E X, u eU and /c, /c* E Nq. 

Note that, because of the particular structure of Fj, with respect to IZ then 
PS is a chain (see [9], page 3) of L differential polynomials, with L as in Section 
3. Let ^ be a characteristic set of [PS] and Aq = An KjX}. By [7] then 

ID = [PS] nK{x} = [A]- 

To compute a characteristic set of [PS] we will use the reduced Groebner basis 
of (PS) with respect to lex monomial order induced by the ranking TZ* . We 
call Q the Groebner basis associated to the system T^(X,U). For that purpose 
we apply the algorithm given in [1], Theorem 6, that we briefly include below 
for completion. 

Given P E K{X U f/}, the lead of P is the highest derivative present 
in P w.r.t. 7^*, we denote it by lead(P). Given P,Q E K{X U U} we de- 
note by prem(P, Q) the pseudo-remainder of P with respect to Q, [9], page 7. 
Given a chain A = {Ai, . . . , At} of elements of ]K{X U U} then prem(P, A) = 
prem(prem(P, At), {Ai, . . . , At-i}) and prem(P, 0) = P. 

Algorithm 4.1. Given the set of differential polynomials PS the next algorithm 
returns a characteristic set A of [PS]. 

1. Compute the reduced Groebner basis Q o/ (PS) with respect to lex mono- 
mial order induced by TZ* . 

2. Assume that the elements of Q are arranged in increasing order Bq < 
Bi <■■■ < Bl_i w.r.t. TZ*. Let A^ = {Bq}. For i from 1 to L - 1 do, if 
lead(Pi) ^lead(Pi_i) then A' := A'-^ U {pTem{Bi, A'~^)}. A := A^~\ 

We are dealing with a linear system of polynomials and computing a Groeb- 
ner basis is equivalent to performing gaussian elimination. Some details on this 
computation were given in [13] and we include them below to be used further 
in this paper. Let M2L be the L x (2L) matrix whose k-th row contains the 
coefficients, as a polynomial in ]K[A:'][V], of the {L — k + l)-th polynomial in 
PS and where the coefficients are written in decreasing order w.r.t. TZ*. 
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ai 



M, 



2L 



ai 



d 



N-o„ 



The polynomials corresponding to the rows of the reduced echelon form 
E2L of M2L are the elements of the Groebner basis associated to V{X, U). Let 
Qo = gn K{X}. From [13], Lemma 20(1): 

Lemma 4.2. The Groebner basis associated to the system V{X, U) is a set of 
linear differential polynomials Q = {Bq, Bi, . . . , Bl^i} where Bq < Bi < ■ ■ ■ < 
Bl_i with respect to the ranking TZ* and Bq & Qq. Hence (PS) fl K{X} = (^0) 
is nonzero. 



4.1 On linear differential polynomials in ID 

The next definitions will play an important role throughout the paper. Given 
a nonzero linear differential polynomial B in [PS] there exist differential oper- 
ators J^i G K[9], i = 1, . . . ,n such that 

n 

B{X,U)=Y,HHX.U)). 

i=l 

If B belongs to ID = [PS] n K{X} then 

n n 

B = 5^ J-.(x, - ai) and J^^iW^)) = ^■ 

i=l 1=1 

If B belongs to (PS) then deg(J^j) < — Oj — 7, i = 1, . . . , n. 

Definition 4.3. Given a nonzero linear differential polynomial B in {PS) with 
B{X, U) = YJUHHX. U)), e K[d]. 

1. We define the co-order of B in (PS) as the highest positive integer c{B) 
such that d^'^^^B e (PS). Observe that 

c{B) = min{A^ — Oj — 7 — ord(i?, Xi) \ i = 1, . . . , n}. 

2. For i E {1, . . . ,n} let ai be the coefficient of d^~"^~^~'^^^^ in Ti. We call 
(«!, . . . , a„) the leading coefficients vector of B in (PS) and we denote it 
by\{B). 
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Let 5* be the leading matrix of the system V{X,U). Denote by S'^ the 
transpose matrix of S. 

Remark 4.4. The i-th row of S equals l{Fi{X,U)), i G {l,...,n}. Given a 
nonzero linear B G {PS), the i-th row of S also consists of the coefficients 
ofUn-,N-^,-^-c(B), J e {!,..., n-1} zn9^— ^-^(^)i^,(X,f/). Therefore if 
oid{B,Uj) <N--fj-j- c{B) then \{B)S = 0, that is 1(5) G Ker(^^). 

Definition 4.5. Given a nonzero linear differential polynomial B in ID with 

1. We define the ID-content of B as the greatest common left divisor of 
J^i, . . . ,Tn (we write gcld(J-'i, . . . , ). We denote it by IDcont(i?). 

2. There exist Ci G K[9] such that J-^ = IDcont(S)£j, i = l,...,n and 
Ci, . . . ,Cn are coprime (we write (£i, . . . , £„) = !)■ We define the ID- 
primitive part of B as 

n 

IDprim(5)(X, U) = -a^). 

i=l 

3. //IDcont(5) G K then we day that B is ID-primitive . 

U A = IDprim(5) = Y^^^^^Ciixi - a^) then c{A) > deg(IDcont(5)) and 
deg(£i) < N -Oi-'j - c{A), i = l,...,n. 

Lemma 4.6. Given a linear differential polynomial B G (PS) fl K{X} then 
IDprim(5) G (PS) nK{X}. 

Proof. For i = 1, . . .n and j = 1, ... n — 1, there exist differential operators 
A, G K[d] such that H,{U) = E;=i' If B{X,U) = EtiH^i ' «0 
then Yl'i=i ^i{.HiiU)) = 0. As a consequence, for each j G {1, . . . , n — 1} then 
Yri=i^ii.'^iji.Uj)) = 0. Let £ = IDcont(5) then = Cd with G K[d] 
and IDprim(i?) = X^iLi ^(a^i — o-i)- Thus C Y17=i ^i^ij = ^-^id £ 7^ so the 
differential operator Yll=i^i^ij — 0- We conclude that E"=i £i(-f^j(f^)) = 0. 
Therefore IDprim(i?) G K{X} which proofs the lemma. □ 

5 Conditions for dim(ID) = n — 1 

Let P{X,U), Fi, Hi be as in Section 2. Let PS = PS(Fi, . . . , F„) and let 
ID be the implicit ideal of P{X,U). Let S and Ml_i be the leading and 
principal matrices of P{X, U) respectively, as defined in Section 3.1. Let Q be 
the Groebner basis associated to the system P{X,U), Qo = Q (1 K{X} and 
denote by \Qo\ the number of elements of Qo- 
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By Lemma 5.3 the ideal (PS) fl K{X} = (Qq) is nonzero. Given a nonzero 
ID-primitive differential polynomial A in (PS) fl K{X} and co-order c{A) as 
defined in Section 4.1, in this section we provide necessary and sufficient condi- 
tions on S, M{L) and A for A{X) = to be the implicit equation of V{X, U). 

5.1 Necessary conditions for dim(ID) = n — 1 

For = ^"=1 Oj = then P(X, U) in a system of n linear equations in n — 1 
indeterminates. 

Lemma 5.1. If N = then dim ID = n — 1 if and only z/rank(5') = n — 1. 

Proof. The matrix M{L) is the n x n matrix whose n x {n — 1) principal 
submatrix is S and whose last column contains Xj — in the i-th row, i = 
1, . . . ,n. Then for linear Ui G K{X} 

rank(5) =n-l^g = {Bq, m - Ui{X), . . . ,Un - Un{X)}. 

Equivalently {-Bo} is a characteristic set of ID. □ 

We prove in the next theorem that one statement is true in general. 

Theorem 5.2. Let V{X, U) he a system of linear DPPEs with implicit ideal 
ID and leading matrix S. //dim ID = n — 1 then rank(S') = n — 1. 

Proof. For = the result was proved in Lemma 5.1. Let us suppose that 
> 1. Let PS = . . . , I z = 1, . . . , n}. Let as defined 

in Section 3 and let X = {xi, . . . , XiN-o,--y-i \ i = 1, . . . n}. Let M2Lh be the 
X (2L^ + 1) matrix whose k-th row contains the coefficients, as a polynomial 

in K[A'][V''], of the {L'^ — l)-th polynomial in PS and where the coefficients 

are written in decreasing order w.r.t. TZ*. Then rank(M2ih) = since M2Lh 

has the L'* x identity matrix submatrix. 

Since rank(5') < n — 1 then Ker(S'-'") 7^ {0}. Given a = (ai,...,a„) G 

Ker(^'^), a ^0 then 

n n 

B{X, U) = Y, a.d^'-^^-'^FiiX, U)-Y, a^d^'-^'-'^ixi - a^) G K[X][V% 

i=l 1=1 

Since rank(M2x,h) = and the first columns of M2Lh correspond to the 
columns indexed by then there exists Qi G K[d] with deg(^j) < A^— Oj— 7— 1, 
i = 1, . . . ,n such that 

n 

B{X,U) +J2S^iF^iX,U)) eK[X]. 

1=1 
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Let 



n 



n 




fN-Oi-j 



1=1 



1=1 



Then G (P5')n]K{X} has leading vector l{Aa) = a and co-order c{Aa) = 0. 

Since dim ID = n — 1 then G ID = [A] for some differential polynomial 
A G which is linear by Lemma 5.3. There exists G K[(9] with 

deg(r'Q,) = c{A) such that Aa = Va{A). Then there is a nonzero G K such 
that l{Aa) = cJ{A). 

Given /3 G Ker(S'^), ^ (3 ^ a there exist G ID = [A] with l{Ai3) = (3, 
c{Ai3) = and a nonzero C/3 G K such that /(A/?) = c^/(v4). Hence /3 = (c_g/ca)a 
which proves that the dimension of Ker(S'"^) is equal to 1. Equivalently, 
rank(5') = n — 1. □ 

As expected, for > the next results show that we need some more 
requirements for dim (ID) = n — 1. 

Lemma 5.3. 1. The number of elements ofQo equals \Qo\ = L—Ta,nk{ML-i)- 

2. For every nonzero linear B G (^o) then \Qq\ > c{B) + 1. 

Proof. 1. Let M2L be the L x 2L matrix defined in Section 4 and E2L its 
redured echelon form. The number of elements of Qq is the number of 
rows in E2L with zeros in the first L — 1 columns. Thus \Qq\ = L — 
rank(ML_i). 

2. Given B G (^^o) = (PS) n K{X}. By definition of c{B) then dB,. . . 
, d^^^B G (PS) nK{X}. Also, there exists k G {!,..., n} such that 
OTd{B,Xk) = N — Ok — '~f — c{B), we can assume that the coefficient of 
XkN-Ok--y in i? is 1. Thus M2L is row equivalent to an L x 2L matrix with 
d^^^^B, . . . , dB, B in the last c{B) + 1 rows. Namely, replace the row of 
M2L corresponding to the coefficients of d^~°^~^~^Fk by d^^^^~^B, t = 
0, . . . , c{B) and next reorder the rows of the obtained matrix. Therefore 
IGol > c(fi) + l. 



If the dimension of ID is n - 1 then ID = [A] for some A G (PS) n K{X} = 
(Qq) and {A} is a characteristic set of ID. By Lemma 5.3 then A is linear and 
we give some more requirements for A in the next theorem. 

Theorem 5.4. Let Q he the Groehner basis associated to the system V{X, U) 
with implicit ideal ID, Qo = Q r\K{X}. //dim ID = n — 1 then ID = [A] where 
A is a nonzero linear differential polynomial verifying. 

1. A is an ID-primitive differential polynomial in (Go). 



□ 
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2. \g^\=c{A) + l. 

Proof. 1. Let C = IDcont(A) and A' = IDprim(/l). If A is not ID-primitive 
then deg(£) > 1 and A = C{A') contradicting that {A} is a characteristic 
set of ID. 

2. If dimID = n-1 and |^o| = m + 1 then Qq = {Bo,Vi{Bo), . . . ,V^{Bo)} 
with Vi e K[d], deg(r'i) > 0, deg{Vi+i) > det(Pi), i = l,...,m. Then 
m < c{Bo). Thus by Lemma 5.3 then \Qq\ = co{Bo) + 1. On the other 
hand ID = [Bq] and A G ID then A = V{Bq) for some V G K[a] 
but we just proved that A is ID-primitive then "D G K. This implies 

\go\ = c{A) + 1. 

□ 

Observe that, if dim ID = n — 1, given A and B nonzero linear ID-primitive 
differential polynomials in (^o) then ID = [Al\ = [B] and c{A) = c{B). 

Corollary 5.5. Let Q he the Groehner basis associated to the system V{X, U) 
with implicit ideal ID, Qq = Q P\ K{X}. //dim ID = n — 1, for every nonzero 
linear ID-primitive differential polynomial A in (Qq) then ID = [A] and \Qq\ = 
c(A) + l. 

5.2 Sufficient conditions for dim(ID) = n — 1 

Given a nonzero linear ID-primitive differential polynomial A in (PS) fl K{X} 
we give sufficient conditions for A to be a characteristic polynomial of ID. 

Lemma 5.6. Let Q he the Groehner hasis associated to the system V{X,U) 
with implicit ideal ID, Qo = Q K{X}. Let S he the principal matrix of 
V{X,U). Given a nonzero linear ID-primitive differential polynomial A in 
[Go) with \Qo\ = c{A) + 1 the following statements hold. 

1. For j = 0, 1, . . . , c{A) there exist Vj G ]K[d] with deg(Pj) = j such that 
Qo = {Bo = Vo{A), 5i = . . . , B,^A) = -D.^aM)}. 

2. //rank(S') = n — 1. Given B G Q\Qo, let us suppose there exists a positive 
integer cb such that 1 < < c{A) and ord(i?, uj) < N — '-/j — — cb, 
j = 1, . . . ,n — 1. Then there exists B G (^o) such that c{B — B) > e^. 

Proof. 1. Since \Qq\ = L — rank(Mi„i) = c{A) + 1, there exists an eche- 
lon form E of M2L whose last c{A) + 1 rows contain the coefficients of 
d^'^^A, ...,dA,A. Then the last c{A) + 1 rows of the reduced echelon 
form E2L of M2L contain the coefficients of -Bc(yi) = 1^c(A){A), . . . ,Bi = 
Vi{A),Bo = Vo{A) for some Vj G K[d] with deg{Vj) = j where j = 
0, 1, ... , c{A). Then Go = {Bq, Si, ... , S,(^)}. 
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2. Let s G {0, . . . , — 1}. Then by 1 the co-order of Bc(a)-s equals 
c{Bc(A)s) = s. Since Bc(a)~s e K{X} then ord(5c(A)-s, Wj) < N - 
Ij — 1 — s, hence by Remark 4.4 then 1(_Bc(a)-s) £ Ker(S'"^). Given 
-B G ^\^o we will prove by induction on s that for s = 0, . . . , — 1 
there exists Cs G (^o) such that c{B — Cg) > s + 1. Then B = Ceg-i- 

There exists Qi G K[d] with deg(^i) < N — Oi — 'j, i = 1, . . . ,n such 
that B = YJ'i=iGi{HX.U)). Let A be the coefficient of d^-°^-^< in 
Qi. By assumption ord(-B,Uj) < N — •jj — •y — cb < — 7j — 7 so 
/3 = . . . , Pn) £ Ker(S'^). Now rank(S') = n — 1 which means that 
dimKer(S'"^) = 1 so there exists G K such that f3 = /i/(i?c(A))- Let 
Cq = fiBc(A) then c{B — Cq) > 1. This proves the claim for s = 0. 

Assuming the claim is true for s — 1, s > 1 then there exists Cs^i G (^0) 
such that c{B - > s. Then B - = J21=i Gt{Fi{X, U)) where 

G K[d] with deg(^f) < iV - o, - 7 - s, z = 1, . . . , n. Let /3f be 
the coefficient of (^^-o»-7-s By assumption ord(i? — Cs-i,Uj) < 

N-y.-y-es < N--fj--f-s so (3' = {l3(, . . . G Ker(5^). Now there 
exists G IK such that /3* = /ij,/(-Bc{A)-s)- Let = Cs-i + iisBc(a)-s 
then c(5 - C,) > s + 1. 

□ 

Theorem 5.7. Let Q he the Groebner basis associated to the system V{X, U) 
with implicit ideal ID, Qq = Q r\ K{X}. Let S be the principal matrix of 
V{X, U). //rank(S') = n—1, for every nonzero linear ID-primitive differential 
polynomial A in (Qo) with \Qo\ = c{A) + l then A is a characteristic polynomial 
of ID. 

Proof. If c{A) = then rank(M2,_i) = L — 1. By Theorem 18(2) and Lemma 
20(4) in [13] and Lemma 5.6(1) then A is a characteristic polynomial of ID. 
Let us suppose that c{A) > 0. For j G {1, . . . , n — 1} let 

Gj = {Beg\ lead(5) = ujk, A; G {0, 1, . . . , iV - 7, - 7}}. 

Then g\go = U^llGj. We use Algorithm 4.1 to prove the result. Set = 
{Bq} then by Lemma 5.6(1) we have prem(i?j, ^*~^) = 0, i = 1, . . . ,c{A), 
thus A' = A°. Given i G {c{A) + 1,...,L - 1} then Bi G Qj^ for some 
ji G {1, ... ,n - 1} and if Gj^ n A'~^ = then pTem{Bi, A'~^) i k{X} and 
A' = A'-^ U {prem{Bi,A'~^)}. If i = c{A) + 1 then Gj^ n A'-^ = 0, we will 
prove by induction on i that for i = c{A) + 2, . . . , L — 1 if fl A^^^ ^ then 
prem(i?j, ^*~^) = and so A^ = A^~^. This proves that ^0 = A^ and A is a 
characteristic polynomial of ID. 

For i = c{A) + 2 let us suppose that Gji H A'^^ ^ then \Q(,\ = L — 
rank(Mi;,_i) = c{A) + 1 implies 

1 < Ci = OTd{Bi,Uj^) - ord(i?i_i,MjJ < c{A). 
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Hence ord(i?j_i, Mj-J < — 7j. — 7 — and by Lemma 5.6(2) there exists 
B'i G (Go) such that c(5i_i - ^) > e^. T\m^d\Bi_i - Bl) e (PS) for 
t = l,...,ej. Therefore = prem(i?j, — Bi) G (PS) with lead(Cj) < 
lead(Ei_i) so Q G (PS) n K{X}. Since prem(5i_i - W^A'-^) = then 
prem(5i,^^-^) = prem(Ci, ^*-i) = prem(Ci,^°) = 0. 

Given i e {c{A) + 3, . . . , L - 1} with Qj^ n A'^'^ 7^ then there exists 
B G {Bc(A)+i, ■ ■ ■ , Bi^i} such that 

1 < Cb = ord(i?j, UjJ - ord(i?, n^J < c{A). 

Hence ord(-B, MjJ < N—'-fj-—^—eB and by Lemma 5.6(2) there exists -B G (^0) 
such that c{B -B)> 63. Thus d^^^B -B)e (PS). Therefore Q = Bi - 
d^'^iB.^i-B) G (PS) with lead(Ci) < lead(Si) then Q = 70^0 + ■ ■ ■+7,_i5,„i 
with 7o, . . . , 7j_i G K. By induction on i then 

prem(Si„i - 5, A'^'^) = prem(prem(5,_i - B, prem(5i„i, A'^^)),A'~^) = 

so we have prem(i?j, ^*~^) = prem(Cj, ^*~^). By induction we also have 
prem(Eo, A''^) = 0, . . . , prem(5i_i, A''^) = then prem(Ci, A''^) = 0. □ 

Corollary 5.8. Given a system V{X, U) of linear DPPEs with implicit ideal 
ID. Let S and Ml_i be the leading and principal matrices ofV{X, U) respec- 
tively. The following statements are equivalent. 

1. The dimension of YD is n — 1. 

2. rank(S') = n — 1 and there exists a nonzero linear ID -primitive differential 
polynomial A such that L — rank(ML_i) = c{A) + 1. 

3. rank(S') = n — 1 and for every nonzero linear ID-primitive differential 
polynomial B then L — rank(ML_i) = c{B) + 1. 

In such situation A{X) = is the implicit equation ofV{X, U). 

Proof. By Theorem 5.2 and Theorem 5.4 then (1)^(2) and by Corollary 5.5 
(1)^(3). Theorem 5.7 proves (3)^(1). □ 



6 Linear perturbations of V{X, U) 

Let V{X, U), Fi, Hi be as in Section 2. Let p be a differential indeterminate 
over K such that d{p) = 0. Denote by Kp = ]K(p) the differential field extension 
of K by p. A linear perturbation of the system V{X, U) is a new system 



V^{X,U) = 




= Pi{U)+pMU) 

= Pn{U)+pMU). 
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where the linear perturbation = (0i(f/), . . . ,(j)n{U)) is a family of hnear dif- 
ferential polynomials in K{f/}. For i = 1, . . . , n let 



Ff{X,U) = F,{X,U)-p4>i{U) and Hf{U) = H,{U) - p <Pi{U) . 

Let PS^ = PS(F/', . . . , F^), a set of linear differential polynomials in [V] C 

Kp{XUf/} and let (PS<^) be the ideal generated by PS^ in Kp[A'][V]. We prove 
next the existence of a linear perturbation cj) such that 9CRes(-Ff , . . . , F^) ^ 0. 



Let us suppose that o„ > o„_i > . . . > Oi to define the perturbation 



where = 1 if Oj 7^ and = if Oj = 0, for z = 1, . . . , tt,. 

Let us suppose that > 1. We denote by M^{L^) the complete differential 
homogeneous resultant matrix for the set of linear differential polynomials 
Hf, . . . , H^. Then M^[L^) is an x matrix with elements in K[p] and 
there exists an x matrix with elements in IK such that M^[L^) = 
M{L^) -pM^. Then 



As expected, the linear perturbation that makes dCIies{Ff , . . . , F^) 7^ is not 
unique. Our proposal can be viewed as a generalization of the characteristic 
polynomial of a matrix to the linear differential case in the spirit of [2] . Namely, 
for On = 1, On-i = ■ ■ ■ = Oi = then we obtain the characteristic polynomial 
of M(L'^), namely 9CRes'^(ifi - p(j)i, . . . , - p(f)n) = det(M(L'^) - p/^h) 
where I^h is the x L'^ identity matrix. 

Let S"^ be the leading matrix of V^{X, U). For i G {1, . . . n} let S'f be the 
(n — 1) X (ri — 1) matrix obtained by removing the i-th row of S"^. 

Proposition 6.1. Given a system V{X, U) of linear DPPEs and the pertur- 
bation defined by (4). The following statements hold. 

1. The determinant of S"^ is nonzero and it has degree n — 1 in p. 

2- If N > 1 then the linear complete homogeneous differential resultant 
9CRes'* (iff , . . . , H^) is a polynomial in lK[p] of degree and not iden- 
tically zero. 

Proof. 1. Observe that S*^ has p's in the main diagonal. 



0=(0i(f/),...,0„(f/))by 




dCRes''{Hf, ...,H^) = det(M<^(L'^)) = det(M(L'^) -pM^). 
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2. We have det{M^{L^)) = p^'^ det{{l /p)M{L^) - M^). If we set y =l/p 
the matrix obtained from yM{L^) — at y = is M^. We will prove 
that det{M^) ^ and therefore the degree of det{M^{L^)) in p is L^. 

In each row and column of there is at least one nonzero entry and 
at most two. Also, if a column has two nonzero entries at least one of 
them is the only nonzero entry of its row. Therefore we can reorganize 
the rows of to get a matrix which has ones in the main diagonal 
and in every row at most one nonzero entry not in the main diagonal. 
Namely the i-th row of is either the only row of with a nonzero 
entry in the i-th column or the row of with its only nonzero entry in 
the i-th column. Thus A^ is a product of elementary matrices and so it 
has nonzero determinant. 

□ 

Theorem 6.2. Given a system V{X, U) of linear DPPEs there exists a lin- 
ear perturbation such that the differential resultant dCKes{Ff , . . . , F^) is a 
nonzero polynomial in K[p]{X} and det(S'^) ^ 0. 

Proof. Let be the perturbation defined by (4) then det(5'^) 7^ by Proposi- 
tion 6.1. If A^ = the result follows from 

n 

aCRes(Ff , ...,Ft) = Y.i-ir^'' det(5f - a,). 

i=l 

If A^ > 1 then dCRes^{Hf, . . . , i/^) ^ by Proposition 6.1. This is equivalent 
by [13], Theorem 18(2) to dGRes{Ff , . . . , F„^) 7^ 0. □ 

If nonzero then (9CRes(F^^, . . . , F^) is a polynomial in p whose coefficients 
are linear differential polynomials in K{X}. We focus our attention next in 
the coefficient of the lowest degree term in p of 9CRes(-F/', . . . , F^). 

Theorem 6.3. Given a system V{X, U) of linear DPPEs let (p be a linear 
perturbation such that dCRes{Ff , . . . , F^) 7^ and det(5'^) 7^ 0. The following 
statements hold. 

1. There exists a linear differential polynomial P in (PS^) fl Kp{X} with 
coefficients in K[p] and content in K such that an ^ ^ where 1{P) = 
(«!, . . . ,a„). 

2. There exists a G M such that 

a„9CRes(Ff , . . . , F^) = (-1)" det{St)dCRes\Hf , . . . , Ht)P{X). (5) 

burthermore — — ^ — - G 
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Proof. 1. By Lemma 5.3 there exists a nonzero linear differential polynomial 
B e (PS^) nKp{X}. There exists c e Kp such that P{X) = cd<^'>B{X) 
has coefficients in K.[p] and content in K. Observe that co(P) = 0. 
Since P G ]Kp{X} then ord(P, uj) < N — 'jj — j and by Remark 4.4 it 
holds 1{B) e Ker((5'^)^). If a„ = then det(S'^) ^ implies = 0, 
i = 1, . . . ,n. This contradicts that c(P) = therefore a„ ^ 0. 

2. Equation (5) follows from [13], Theorem 18(1). Since the content of 
P{X) belongs to K then a„ divides det(^^)aCRes''(i/f , . . . , H^) in K[p]. 

□ 

Remark 6.4. In the situation of the previous theorem and with the perturba- 
tion (f) defined by (4) we can add the following remarks. For i = 1, . . . ,n by 
formula (3) then det{Sf)dCRes'^{Hf, . . . ,H^) is the coefficient of XiN-Oi-'y in 
9CRes(F/, . . . , F^). By 1 and 2 in the previous theorem then 

an det{St)dCRes^{Ht, . . . , i/^) = det{St)dGRe?!' {Hf , Ht)ai, 

hence det(5f) = det(5^)ai/a„ G K[p]. Let a = gcd(a„, det(5^)). If a E K 
then the degree of det(S'f) in p is greater or equal than n — 2 which is not 
possible then gcd(a„, det(S'^)) G K[p]\K. 

Let Dfi, be the lowest degree of p in (9CRes(F/, . . . , F^) and let A^,^ be the 
coefficient of p^'p in 9CRes(F/, . . . , F^). We call the degree of the perturbed 
system V^{X,U). Observe that 

D^ = 0^ 9CRes(Fi, . . . , F„) ^ 0. 

We write = -1 if dCRes{Ff, . . . , F^) = and so 

D^>0^ dCRes{Ff, ...^F^)^ 0. 

In the remaining parts of this section we assume that 9CRes(F/', . . . , F^) ^ 
0. Let PS = PS(Fi,...,F„) and let ID be the implicit ideal of V{X,U). 
We will use dCRes{Ff , . . . , F^) to provide a nonzero ID-primitive differential 
polynomial in (PS) n K{X}. 

Lemma 6.5. The linear differential polynomial A^^ belongs to (PS) r\K{X} . 

Proof By [13], Proposition 16 then 9CRes(F/, . . . , F^) G (PS<^) n Kp{X} 
and it equals p^'I'^Ad^ + pA') for some A' G KpjX}. Therefore the linear 
polynomial Ajj^ + pA' G (PSp) fl IKpjX} and there exist J-^ G Kp[d] with 
deg(J-i) < AT - - 7 such that Ad,{X) + pA'{X) = Eti H^^f {X , U)) . 
Then AD,iX)+pA'{X) = Y.tl^^{^^ ' and T.ti ^^{Hf {U)) = 0. 
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For each i e {1, . . . ,n} there exists £j G K[d] and J^- G ]Kp[c}] such that 
= A+pJ-/. Then ADj,X)+pA'iX) = Eti ^i) +PEtlJ'^i^^- ^^i) 

and hence Ad^{X) = Yl^=i ^ii^i ~ ^.i)- the other hand 

n n n n 

1=1 «=1 i=l i=l 

which imphes Er=i A(^«(f/)) = 0. Thus Ad^{X) = ^^'L^ f/)) with 

deg(£i) <N-Oi--f then G (PS) n K{X}. □ 

Let A^p be the ID-primitive part of Ad^. We call A^ the differential polyno- 
mial associated to V^{X,U). We relate with c(A0) and give conditions on 
for ^^(X) = to be the imphcit equation of V{X, U). 

Remark 6.6. LetV^{X,U) andV^{X,U) be two different linear perturbations 
ofV{X, U) with degrees > and > 0. Let A^ and A^ be the associated 
differential polynomials. 

1. As illustrated in Example 1 of Section 7 the degrees and may be 
different. 

2. //dim ID = n — 1 then A^ = ^A^ for some 7 G K. 

Theorem 6.7. Let V^{X,U) be a perturbed system of the system V{X,U) 
of degree > 0. Let Q be the Groebner basis associated to V{X,U) and 
go = Qn K{X}. Then |^o| - I < D^. 

Proof. Let M^{L) be the differential resultant matrix of F^, . . . , and Mf_^ 
the L X (L — 1) principal submatrix of M^{L). By equation (2) then 

dCRes{Ft,...,Ft) = Y, E hkdei{MtJ{x,,-d^a,), 

with M^. an (L — 1) x (L — 1) submatrix of M^_-^ and bi^ = ±1 according to the 
row index of Xik — d^ai in the matrix M^{L). Let Mx-,. be the (L — 1) x (L — 1) 
submatrix of the principal matrix Ml^i of V{X,U). Let = rank(M2;.j.), 
then there exists an invertible matrix Ei^ of order L — 1 and entries in K such 
that the last L — 1 — rows of Ei^M^.^ are zero. If we divide each one of the 
last L — 1 — Tjfc rows of Ei^M^^^ by p we obtain a matrix Nf^ such that 

n N—Oi—^ 

det M^{L) = J2 Yl det{Ni){x,k - d'a,). 

i=l k=0 
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This proves that L—l—rik < for alH = 1, . . . , n and k = 0,1, . . . , A^— Oj— 7. 
Now rank(Mi„i) > so 

l^ol -l = L- rank(Mi_i) - 1 < L - 1 - < D^. 

□ 

We showed that > \Qo\ — 1 = L — rank(ML_i) — 1 > c{A^) and in 
general equahty does not hold ( see examples in Section 7). 

Corollary 6.8. Let V{X, U) be a system of linear DPPEs with implicit ideal 
ID and leading matrix S. Let V^{X,U) be a perturbed system of V{X,U) of 
degree > 0. Let be the differential polynomial associated to V^{X, U). If 
rank(S') = n — 1 and = c{A^) then ID has dimension n — 1 and y4^(X) = 
is the implicit equation ofV{X, U). 

Proof. If = c{Afp) then by Theorem 6.7 we have \Qq\ < c{A^) + 1. By 
Lemma 5.3 then \Qq\ = c(A(^) + 1 and by Corollary 5.8 the result follows. □ 

7 Implicitization algorithm for linear DPPEs 
and examples 

Let V{X, U) be a system of linear DPPEs with implicit ideal ID. In this 
section, we give an algorithm that decides whether the dimension of ID is 
n—1 and in the affirmative case returns the implicit equation of V{X, U). For 
this purpose let S and Mi_i be the leading and principal matrices of V{X, U) 
respectively. Let P</,(X, U) be a perturbed system of V{X, U) of degree > 0. 
Let Ad^ be the coefficient of p^"* in dGKe?,{Ff , . . . , F^) and A^f, the differential 
polynomial associated to V^{X, U). 

Algorithm 7.1. • Given the system V{X, U) of linear DPPEs. 

• Decide whether the dimension is n — 1 and in the affirmative case 

• Return a characteristic polynomial of YD. 

1. Compute I ai\k{S). 

2. //rank(5) <n-l RETURN "dimension less than n - 1". 

3. Compute V^{X,U) with perturbation given by (4). 

4. Compute dCRes{Ff, F^), and A^^. 

5. IfD^ = RETURN Aj,,. 
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6. Compute and c{A^). 

7. IfD^ = c{A^) RETURN A^. 

8. Compute Ta,nk{ML_i) . 

9. If L- rank(ML_i) > c{A^) + 1 RETURN "dimension less than n - 1". 
10. IfL- rank(ML_i) = c{A^) + 1 RETURN A^. 

The next examples were computed in Maple. The computation of differ- 
ential resultants was carried out with our Maple implementation of the linear 
complete differential resultant, available at [11]. 



7.1 Example 1 

Let K = Q, 9 = ^ and consider the system V{X, U) of linear DPPEs providing 
the set of differential polynomials in K{xi, X2, xajjMi, M2}, 

Fi{X, U) = Xi + Ui - U2 + Mil - U12 - 4^21 - 3M22, 

F2{X, U) = X2 + U2 + uii - U22, 
Fs{X, U) =X3 + U2 + nil + ^21. 

The set PS(Fi,F2,F3) contains L = 13 differential polynomials and 7 = 0. 
The leading matrix S of V{X, U) has rank 2 and equals 



S 



-3 1 
-1 
1 1 



We consider the perturbation = {(f)i{U) , (f>2{U) , (ps^U)) with 




There exists a differential polynomial P{X) G (PS)nK{X} with coefficients 
in and content in K such that the determinant of the 13x 13 matrix M^[13) 
equals 

dCRes{Ft,Fi,F^) = dCRes'' {Hf , Hi) P{X) = 

p (1 + 4p + - / + 2/ + 11 p3 + / - 12p2 _ 4p6 ^ p^)P{X). 

Then = 1 and the coefficient of p in dCKes{Ff , F2, F^) is 

Ad. = X12 - X2 - 2 X21 - 2 X22 + a;33 + X32 + 0:31 + 0:3. 
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We have A^^ = + C2{x2) + ^^{xs) with 

£2 = -1 - 3 9 - 4 92 - 2 9=^ = -(9 + 1) (2 92 + 2 9 + 1), 
£3 = 1 + 29 + 292 + 29=^ + 9' = (92 + 1) (9 +1)2. 

Therefore £ = gcld(£i, £2, £3) = 1 + d and A^f, = X12 — X2 — 2x2i — 2x22 + 
3^33 + X32 + 2:31 + X3 with c{A^) = 1. Then = c{A^). We conclude that the 
dimension of ID is n — 1 = 2 and its imphcit equation A^{X) = 0. 
If we consider the perturbation = {ipi{U) , ip2{U) , xl^s^U)) with 

{U22 + U1, i = l, 
U2, 1 = 2, 

Mil, ^ = 3. 

There exists a differential polynomial P{X) G (PS)nK{X} with coefficients in 
]K[p] and content in IK such that the determinant of the 13 x 13 matrix M^j(13) 
equals 

9CRes(Fi'^, Ft, F^) = dCRes'^iNf, Ht, Hf)P{X) = 
- p2 (l+p) _ ^-^^5 _ 53^4 _ gQp3 _ 33^2 _ _ i)p(x). 

Then = 2 but the coefficient of ^2 is = Ad^. Therefore A^ = A^ as 
expected. 



7.2 Example 2 

Let K = Q, 9 = ^ and consider the system V{X, U) of linear DPPEs providing 
the set of differential polynomials in Kjxi, X2, 3:3, Xi}{ui, U2, U3}, 

Fi{X, U) = xi - 2ui + U'i- 3u2i + U31, 
F2(X, U) = X2 + 2ui -u-i- uii + 3m22 - M32, 
F3(X, f/) = X3 + 2mi - M3 + 2m2i + M32, 
Fa{X, U) = Xi + 2ui - U3 + 3^21 - 2^31. 

The set PS(Fi, F2, F3, F4) contains L = 18 differential polynomials and 7 = 
7i = 1. The leading matrix 5* of 'P{X, U) has rank 3 and equals 

1 -3 -2 ■ 

-1 3 -1 

10 0' 

-2 3 2 
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We consider the perturbation 



(0l(f/), 02(^7), 03(^7)) With 



I 



U21 + Us, i 

Uu + U2, i 

ui, i 

U31, i 



1, 
2, 
3, 
4. 



There exists a differential polynomial P{X) G (PS)n]K{X} with coefficients in 
K[p] and content in K such that the determinant of the 18 x 18 matrix M<^(18) 



aCRes(Ff , F^, Fi, ) = dCRes^Hf, Ht)P{X) = 

(-2/ - 3944p^ + 78V - 481p^ - 379p^ + 108 + 4484p^ + 212p^ + p^^ 
- 642p + 527p^ - 7p^°)P(X). 

Then = 3 and the coefficient of p^ in aCRes(Ff , ^2^, F^, Ff) is 

^D, = 

972x12 - 864x13 - 972xi4 - 216x22 + 648x32 - 972x33 + 540x42 - 972x44- 

We have Ajj = £i(xi) + £2(2^2) + '^si^s) + £4(x4) with 

£1 = -8649^ + 972^2 - 9729^ = -1089^(85 - 9 + 9(9^), 
£2 = -216a^ 

£3 = -972(9^ + 648(9^ = -32492(39 - 2), 
£4 = 5409^ - 9729^ = -10892(-5 + 99^). 

Therefore £ = —1089^ and = 8xii + 9xi2—9xi+2x2 — 6x3 + 9x31+9x42 — 5x4 
with co(A^) = 2. Then > c{A^). 

Replace p in M^{18) by zero to obtain M(18) whose principal 18 x 17 
submatrix is Ml^i. Compute L — rank(Mi;,_i) = 3. Then c{A^) + 1 = L — 
rank(Mi,_i). We conclude that the dimension of ID is n — 1 = 3 and its implicit 
equation A^{X) = 0. 

7.3 Example 3 

Let K = Q(t), 9 = f and consider the system V{X,U) of linear DPPEs 
providing the set of differential polynomials in K{xi, X2, X3}{-Ui, M2}, 



equals 



-Fi = xi - 3 + Mil + Ui2 - M2 - 4 M21 - 3 M22, 

F2 = X2 + Uii + U2- U22, 

Fs = X3 + 2 + Uu + tu2 + U21. 
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Then the set PS(Fi, F2, -F3) contains L = 13 differential polynomials and 7 = 0. 
The leading matrix S of V{X, U) has rank 2 and equals 



S 



3 1 
1 
1 1 



We consider the perturbation = {(f)i{U) , (f>2{U) , (ps^U)) with 




Ui2 + U2, i 



I = 



I = 



1, 

2, 
3. 



There exists a differential polynomial P{X) G (PS) fl K{X} with coeffi- 
cients in ]K[p] and content in K such that the determinant of the 13 x 13 
matrix M0(13) equals dCRes{Ff , F2 , F^) = pP{X). In this case a^p = 
det{Sf)dCRes'^{Hf,H2,H^) where 1{P) = {ai, 02,(^3) is the leading vector 
of P. Then = 1 and the coefficient Ae>^ of p in dCIies{Ff , F2, F^) equals 



£1 = -64 - 2656 + 912 1 - 8 + (468 + 57 + 33 - - 13 - 522 1) d 

+ (-364 + 79 - 6 1 - 10 - t^) (9^ 

+ (68 1 + 37 -Ut^-t^- 296) a^ 
£2 = (730 1 + 860 -37^3 - 229 + 15t^ + t^)d 

+ (252 1 + 19 + + 1752 - 421 + 17 1^) (9^ 

+ (56 - 272 1 - 148 + 4 + 1184) 9^ 
£3 = - 16 ^2 + 228 1 - 2 - 664 + (-2 + 212 - 436 1 - 18 - 664) 9 

+ (-1856 - 64 + 309 t^-5t^ + 276 1) d'^ 

+ (-2 + 210 1 + 32 - 524 - 32 1^) 

+ (-68 1 - 37 + 14 + + 296) 9^ 

Using the Maple package Ore Tools we check that (£1, C2, ^C^) = 1, then A^^ = 
A^ and co(y40) = 0. Replace p in M^{13) to obtain M (13) whose principal 
13 X 12 submatrix is Ml_i. Compute L— rank(M2,_i) = 2 which is greater than 
c{A^) + 1 = 1. We conclude that the dimension of ID smaller than n — 1 = 2. 
If we take a different perturbation ip = {ipi{U) , ip2{U) , ipsiU)) with 



£i(xi-3) + £2(a;2) + /:3(a;3 + 2) 



with 




U2, i = 2, 

Mil , i = 3. 
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We obtain = 1 and the coefficient of p in IS 

Ad^ = /Ci(Xi - 3) + /C2(X2) + /C3(X3 + 2) 

with 

/Ci = (-76^3 + 154^2 + I2t- 156 + + 

+ (204 1 - 100 + + 8 - 94 1^) 92 

+ (-58 + 12 + 88 1 - 12 + t^) 9^ 
/C2 = -228 1 + 16 + 664 + 2 + (-^5 _ 904 1 + 820 + 62 + 90 - 13 1^) d 

+ (-t^ - 17 + 346 - 892 1 + 412 + 14 1^) 

+ (232 - 352 1 - 4 - 48 + 48) d^, 
/C3 = -228t+ 16^2 + 664 + 2^3 + (Mt^ + 892t - 664 + 2t^ - 244 1^)^ 

+ (-156 -406t2 + 54t^ + 676t + 5t^) 9^ 

+ (-80 + 28 + 2 + 64 + 60 1) 

+ (58 - 88 1 - - 12 + 12) d\ 

Using the Maple package OreTools we check that (/Ci, /C27 f^s) = 1, then Ajj^ = 
A^ and c{A^) = 0. Observe that there is no 7 e IK such that A^f, = 'jA^. 

Acknowledgements. The author thanks J. Rafael Sendra for kindly reading 
and commenting on some parts of this work. 
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